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1) INTRODUCTION: 

As a matter of fact, the solution space of many systems of ordinary differential (OD) or partial 
differential (PD) equations in engineering or mathematical physics " can/ cannot" be parametrized 
by a certain number of arbitrary functions behaving like "potentials" . In view of the explicit ex- 
amples to be met later on, it must be noticed that the parametrizing operator, though often of the 
first order, may be, on the contrary, of arbitrarily large order. Among the well known examples, we 
recall that a classical OD control system is parametrizable if and only if it is controllable (Kalman 
test of 1969 in [9]). Among PD systems, the electromagnetic (EM) field, solution of the first set of 
4 Maxwell equations, admits a well known first order parametrization by means of the ^-potential 
while the EM induction, solution of the second set of 4 Maxwell equations (in vacuum), also admits 
a first order parametrization by means of the so-called A-pseudopotential. On the contrary, it is 
now known that, contrary to the EM situation, the set of 10 second order linearized Einstein equa- 
tions (in vacuum) cannot be parametrized and cannot therefore be considered as field equations 
(see [18,24] for more details; see also [30] and |http://wwwb. math.rwth-aachen.de/O reModules for 
a computer algebra solution). One among the best interesting and useful cases is concerned with 
continuum mechanics where the first order stress equations (in vacuum) admits a rather simple 
second order parametrization by means of the single Airy function in dimension 2 and it is not so 
well known (!) that a much more complicated second order parametrization can be achieved in 
dimension n > 2 by means of n 2 (n 2 — 1)/12 arbitrary functions. 

It is now known that all the above problems are only very particular cases of a much more 
sophisticated (!) and general situation involving the formal theory of systems of PD equations 
pioneered by D.C. Spencer in 1970 [26] (jet theory, diagram chasing, differential sequences, ...)(see 
[22,23] for more details) and differential modules in the framework of " algebraic analysis" pioneered 
about at the same time by V.P. Palamodov [17] for the constant coefficients case, M. Kashiwara [10] 
and B. Malgrange [15] for the variable coefficients case (see [23,24] for more details and also consult 
Zentralblatt Zbl 1079.93001 for a review). The corresponding duality theory, that is at the heart of 
all the previous examples and will be a central tool in this paper, highly depends on (hard) homo- 
logical algebra techniques (localization, resolutions, extension modules,...) which cannot be avoided. 

The purpose of this paper is to apply these techniques in a way as simple and self-contained as 
possible in order to give a positive and explicit answer concerning the possibility to exhibit a first 
order parametrization of the stress/couple-stress equations met in the study of Cosserat media. At 
the same time, as a corollary of the homological test, we shall give for the first time the reason for 
which the compatibility conditions (CC) for the deformation tensor in classical elasticity theory 
are second order while the corresponding CC for Cosserat fields [11] are only first order and explain 
why this order is equal to the order of the corresponding parametrization. 

At the end of the paper, we shall give hints in order to explain why, though the "fields" and 
their CC in classical and Cosserat elasticity theories look like completely different at first sight, 
therefore providing different presentations of the corresponding field equations, nevertheless the 
possibility to obtain a parametrization in one framework necessarily implies the possibility to have 
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a parametrization in the other framework and vice-versa. Though striking it may look like in such 
an engineering background, this totally not evident result, which was not known up to now, is one 
of the simplest consequences of a delicate theorem of homological algebra (see [25], Th 7.22, p 212). 
In particular the reader must look at the first section below with care, even though it does not seem 
to have anything to do with Cosserat media. 

2) MOTIVATIONS: 

In the middle of the last century, commutative algebra, namely the study of modules over 
rings, was facing a very subtle problem, the resolution of which led to the modern but difficult 
homological algebra with sequences and diagrams. Roughly, one can say that the problem was 
essentially to study properties of finitely generated modules not depending on the "presentation" 
of these modules by means of generators and relations. This very hard step" is based on homo- 
logical/cohomological methods like the so-called "extension" modules which cannot therefore be 
avoided [1,7,16,25]. 

EXAMPLE 2.1: In order to sketch this problem, let us present two simple examples. In the 
first case with standard notations, everybody will understand at once that integrating the second 
order OD equation y = with m = n = 1 is equivalent to integrating the system of two first 
order OD equations y 1 — y 2 = 0, y 2 = 0. However, even with m = n = 2 and the same two 
unknowns u, v in both cases, it is not evident at all that integrating the second order PD equation 
di2U — d22V — u = is equivalent to integrating the system of two fourth order PD equations 
dii22U — d 12 22V - d 2 2V - u = 0, d 1112 u - d 112 2V - d\\u = (exercise !). 

As before, using now rings of " differential operators" instead of polynomial rings led to dif- 
ferential modules and to the challenge of adding the word " differ entiaT in front of concepts of 
commutative algebra. Accordingly, not only one needs properties not depending on the presenta- 
tion as we just explained but also properties not depending on the coordinate system as it becomes 
clear from any application to mathematical or engineering physics where tensors and exterior forms 
are always to be met like in the space-time formulation of electromagnetism. Unhappily, no one of 
the previous techniques for OD or PD equations could work !. 

By chance, the intrinsic study of systems of OD or PD equations has been pioneered in a totally 
independent way by D. C. Spencer and collaborators after 1960 [26], using jet theory and diagram 
chasing in order to relate differential properties of the equations to algebraic properties of their 
symbol, a technique superseding the " leading term" approach of Janet in 1920 [8] or Grobner in 
1940 [6] but quite poorly known by the mathematical community, even today. 

Accordingly, it was another challenge to unify the "purely differential' approach of Spencer 
with the "purely algebraic" approach of commutative algebra, having in mind the necessity to use 
the previous homological algebraic results in this new framework. This sophisticated mixture of 
differential geometry and homological algebra, now called " algebraic analysis" , has been achieved 
after 1970 by V. P. Palamodov for the constant coefficient case [17], then by M. Kashiwara [10] 
and B. Malgrange [15] for the variable coefficient case. 

Let k be field of characteristic zero and \ = (xi, Xn) be indeterminates over k. We introduce 
the ring A = k[xi, ■■■,Xn] of polynomials with coefficients in k. When a given system of linear PD 
equations of order q is given, it defines by residue a differential module M over the underlying 
ring D of differential operators. More precisely, let us introduce n commuting derivatives d\,...,d n 
for which k should be a field of constants and define the ring D = k[d] — k[di,...,d n ] of differ- 
ential operators with coefficients in k. Then D and A are isomorphic by di <-> Xi- However, the 
(non-commutative) situation for a differential field K with subficld of constants k escapes from the 
previous (commutative) approach and must be treated "by its own" . For this, let \i = fjiii, ...,/U„) 
be a multi-index with length | /i |= + We set fi+li = (fi±, /Uj-i, Atj + 1, Mi+ij ■••> Mn) and 

we say that \x is of class i if n\ = ... = Hi-i = 0,/ij ^ 0. Accordingly, any operator P = a^d^ € D 
acts on the unknowns y k for k = 1, m as we may set d^y k = y k with yft = y k and introduce the 
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jet coordinates y q = {y^ \ k = 1, ...,m;0 <| < q}. It follows that, if a system of PD equations 
can be written in the form $ T = a^y^ = with a £ if, we may introduce the filtrcd differential 
module M = Dy/D<&. Then we define the (formal) prolongation of <I> T with respect to di to be 
di$ T = a^ M y^ +1 . + dia^yk and induce maps di : M — > M : y^ — > by residue. 

First of all, we sketch the technique of " localization" in the case of OD equations, comparing 
to the situation met in classical control theory where n = 1. For this, setting as usual d = d\ = 
d/dt = dot, we may introduce (formal) unknowns y 1 , ...,y m and set Dy = Dy 1 + ... + Dy m ~ D m . 
If we have a given system $ = of OD equations of order q, a basic question in control theory is 
to decide whether the control system is " controllable" or not. It is not our purpose to discuss here 
about such a question (see [23,24] for more details) but we just want to state the final formal test 
in terms of a property of the differential module M = Dy/D$. Care must be taken that in the 
sequel, for simplicity and unless needed, we shall not always put a " bar''' on the residual image of 
y in the canonical projection Dy — ► M. We explain our goal on an example. 

EXAMPLE 2.2: With m — 3 and a constant parameter a, we consider the first order system 
$ x = y 1 — ay 2 — y 3 = 0, $ 2 = y 1 — y 2 + y 3 = 0. Any engineer should want to apply Laplace 
transform y(s) = e st y(t)dt to this system. However, using the integration by part formula 
'e st y{t)dt — [e s *y(t)]fj° — sy(s) we should eventually need to know y(0) though the Kalman test 
of controllability is purely formal as it only deals with ranks of matrices [9]. Since a long time 
we had in mind that setting y(0) — was not the right way and that Laplace transform could be 
superseded by another purely formal technique. For this, let us replace "formally" d by the purely 
algebraic symbol \ whenever it appears and obtain the system of linear equations : 

l 2 3 n 1 1 2 3 1 x(x + a) 3 2 X(X + 1) 3 

Xy - ay z - XV = 0, V - XV + XV = => y L = V , V = ~\ V 

X -a x - a 

but we could have adopted a different choice for the only arbitrary unknown. At this step there 
are only two possibilities : 

•a^0,l =>no "simplification" may occur and, getting rid of the common denominator, we get an 
algebraic parametrization leading to a differential parametrization as follows: 

V 1 = X(X + a)z, y 2 = x{x + IK y 3 = (x 2 - a)z y 1 = z + az, y 2 = z + i, y 3 = z - az 

•a = 0ora=l=>a "simplification" may occur and no parametrization can be found. For exam- 
ple, with a = 0, setting z = y 1 — y 3 we get x z — that is to say z = 0. 

Recapitulating, we discover that a control system is controllable if and only if one cannot get 
any autonomous element satisfying an OD equation by itself. For understanding such a result in 
an algebraic manner, let M be a module over an integral domain A containing 1. A subset S d A 
is called a multiplicative subset if 1 £ S and Vs,i £ S => st £ S. Moreover, we shall need/use the 
Ore condition on S and A, namely aS PI sA ^ 0, Va £ A, s £ S. 

DEFINITION 2.3: S^A = {s'^s £ S,a £ A/ ~} with s^a - t^b <=> 3u, v £ A, us = vt £ 
5, ua = vb. 

Next, for any module M over A, we dehne S _1 M = S^A® A M and t 5 (M) = {x £ M \ 3s £ 
S, sx = 0} in the exact sequence — > ts(M) — > M — > S' _1 M where the last morphism is x — > l _1 a;. 

EXAMPLE 2.4: 5 = A - {0} => S"" 1 ^ = Q(A) field of fractions of A and we introduce the 
torsion submodule t s (M) = t(M) = {x £ M | 30 7^ a £ A, ax = 0} of M. 

In the case of a torsion-free module, reducing to the same denominator as in the control exam- 
ple, we obtain the following classical proposition amounting to exhibit a parametrization [23,25]. 
However, the reader must notice that it is unuscfull in actual practice as one needs a test (like the 
Kalman test) for checking the torsion- free condition. This will be the hard part of the job !. 
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PROPOSITION 2.5: When M is finitely generated and t(M) = 0, from the inclusion M C 
Q(A)®aM , we deduce that there exists a finitely generated free module F with M C F. 

EXAMPLE 2.6: In 2-dimensional classical elasticity, let us consider the well known stress equa- 
tions: 

8 1 a ll +d 2 a 21 =0,d 1 a 12 +d 2 <J 22 =0 

where we must keep in mind that the stress tensor density is symmetric, that is a 12 — a 21 . Re- 
placing di by Xii we mav localize and obtain: 

Xi^ U + X2<x 21 = 0, xi^ 12 + X2<r 22 = 
Reducing the fractions to the same denominator, we get: 

11 X2 2i (X2) 2 12 22 Xi 12 
a = a = a , a = a 

Xi X1X2 X2 

and obtain therefore the subvector space over Q(xi, X2) : 

a 11 = (X2) 2 0,<T 12 = a 21 = -XiX20^ 22 = (Xift 
a result providing at once the well known parametrization by the Airy function: 

C 11 = <9 2 2</>, = 0-21 = ^^120, CT22 = <9ll</> 

It may be interesting to compare this purely formal approach to the standard analytic aproach 
presented in any textbook along the following way. From the first stress equation and Stokes 
identity for the curl, there exists a function ip such that a 11 = c^, cr 21 = —d\tp. Similarly, from 
the second stress equation, there exists a function ip such that a 22 — d\ip, a 12 = —82^- Finally, 
from the symmetry of the stress, there exists a function 4> such that ip — d2(f>,ip — di4> an d we 
find back the same parametrization of course. The reader must notice that, in this example, one 
can check that the parametrization does work but no geometric inside can be achieved in arbitrary 
dimension n > 2 even though exactly the same procedure can be applied. 

EXERCISE 2.7: Apply similarly the localization technique in the case of the two sets of Maxwell 
equations and compare to the standard analytic approach to be found in any textbook. 

The extension of the above results to the non-commutative case D — K[d] where if is a differ- 
ential field with n commuting derivations di,...,d n can be achieved but is more delicate [ 23]. 

EXERCISE 2.8: When a — a(t) in Example 1.6, we let the reader prove that the controllability 
condition is now the Ricatti inequality a + a 2 — a^Oina coherent way with the constant coefficient 
case already considered. 

Taking into account the works of Janet and Spencer, the study of systems of PD equations cannot 
be achieved without understanding involution and we now explain this concept. For this, changing 
linearly the derivations if necessary, we may successively solve the maximum number of equations 
with respect to the jets of order q and class n, class (n — 1),..., class 1. At each order, a certain 
number of jets called principal can therefore be expressed by means of the other jets called para- 
metric. Moreover, for each equation of order q and class i, d\,...,di are called multiplicative while 
di + \,...,d n are called nonmultiplicative and d\,...,d n are nonmultiplicativc for all the remaining 
equations of order < q — 1 . 

DEFINITION 2.9: The system is said to be involutive if each prolongation with respect to a 
nonmultiplicative derivation is a linear combination of prolongations with respect to the multi- 
plicative ones. Using Spencer cohomology, one can prove that such a definition is in fact intrinsic 
[22,23,26]. 



(Xi) 2 
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EXAMPLE 2.10: The system j/n — 0, 2/13 — y<i = is not involutive. Effecting the permutation 
(1, 2, 3) — *■ (3, 2, 1), we get the new system 2/33 = 0, 2/13 - yi = 0. As dij/ 33 - d 3 (yi 3 - 2/2) = 2/23 and 
<ii2/23 - ^2(2/13 - 2/2) = 2/22, the new system 2/33 = 0, 2/23 = 0, 2/22 = 0, 2/13 - 2/2 = is involutive with 
1 equation of class 3, 2 equations of class 2 and 1 equation of class 1. 

EXAMPLE 2.11: The Killing system yf + y l j = is not involutive but the first prolongation 
y\ + y l j = 0, y^j — is involutive. This is the reason for which the Ricmann tensor is a first order 
expression in the metric and Christoffcl symbols and thus second order in the metric alone (for 
more details, see [22], p 249-258). 

APPLICATION 2.12: t(M) = M if and only if the number of equations of class n is to. Oth- 
erwise there is a strict inclusion t(M) c M. 

PROPOSITION 2.13: ([22,26]) The following recipe (already used in the example at the begin- 
ning of this section) will allow to bring an involutive system of order q to an equivalent (isomorphic 
modules) involutive system of order 1 with no zero order equations called Spencer form: 

1) Use all parametric jets up to order q as new unknowns. 

2) Make one prolongation. 

3) Substitute the new unknowns. 

3) GROUP FOUNDATION OF MECHANICS: 

This section, which is a summary of results already obtained in [21,22], is provided for fixing the 
notations and the techniques leading to various (linearized) differential sequences. All the results 
presented are local ones. A corresponding non-linear does exist but is out of the scope of this paper 
[12,21]. 

Let A be a manifold of dimension n with local coordinates x — (x 1 ,...,x n ) and latin indices 
i,j = 1, n. We denote by T = T(X) the tangent bundle to A and by T* = T*(A) the cotangent 
bundle to A while A r T* is the bundle of r- forms on A. Also, we denote by J q (T) the q-jet bundle 
of T, that is to say the vector bundle over A having the same transition rules as a vector field and 
its derivatives up to order q under any change of local coordinates on A. Let now G be a Lie group 
of dimension p with identity e, local coordinates a — (a , a p ) and greek index r. We denote by 
Q = T e (G) the corresponding Lie algebra with vectors denoted by the greek letters A. As usual, 
we shall identify a map a : X — > G called gauging of G over A, with its graph I-tlxG which 
is a section of a trivial principal bundle, and, similarly, use the same notation for a bundle and its 
sheaf of (local) sections as the background will always tell the right choice. In particular, when 
differential operators are involved, the sectional point of view must automatically be used. Such a 
convention allows to greatly simplify the notations at the expense of a slight abuse of language. 

DEFINITION 3.1: A Lie group of transformations of a manifold A is a lie group G with an 
action of G on A better defined by its graph A x G — > A x A : (x, a) — > (x, y = ax = f(x,a)) 
with the properties that a(bx) = (ab)x and ex = x, Vx e A, Va, b 6 G. 

It is sometimes useful to distinguish the source x from the target y by introducing a copy Y of 
A with local coordinates y = (y 1 , ...,y n ). Such groups of transformations have first been studied 
by S. Lie in 1880. Among basic examples when n = 1 we may quote the affine group y = ax + b 
and the projective group y = (ax + b)/(cx + d) of transformations of the real line. When n = 3 we 
may quote the group of rigid motions y — ax + b where now a is an orthogonal 3x3 matrix and b 
is a vector. Such a group is known to preserve the euclidean metric lo = (lu ij = and thus the 
quadratic form ds 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 — LO % ^dx % dxK When n = 4 we may quote the con- 
formal group of space-time with 15 parameters (4 translations, 6 rotations, 1 dilatation, 4 elations) 
preserving the Minkowski metric ui or the quadratic form ds 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 — c 2 (dt) 2 
up to a function factor, where now c is the speed of light and t the time. Among the subgroups, 
we may consider the Weyl group with 11 parameters preserving ui up to a constant factor and the 
Poincare group with 10 parameters preserving w. 
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Only ten years later, in 1890, S. Lie discovered that the Lie groups of transformations were 
only examples of a wider class of groups of transformations, first called infinite groups but now 
called Lie pseudogroups. 

DEFINITION 3.2: A Lie pseudogroup T of transformations of a manifold X is a group of 
transformations y — f(x) solutions of a (in general nonlinear) system of PD equations, also called 
system of finite Lie equations. 

Setting now y = x + t£(x) + ... and passing to the limit for t -> 0, that is to say linearizing the 
defining system of finite Lie equations around the q-jet of the identity y = x, we get a linear sys- 
tem R q C J q {T) for vector fields, also called system of infinitesemal Lie equations, with solutions 
0CT satisfying [6, 0] C 6. It can be proved, for the same testing type reasons, that such a sys- 
tem may be endowed with a Lie algebra bracket on sections £ q : (x) — > (x, £ fe (x), £f(x), £^-(x), ...) 
that we shall quickly define (see [22] for more details and compare to [12]). Such a bracket 
on sections transforms R q into a Lie algebroid and we have [R q ,R q ] C R q . Let us first de- 
fine by bilinearity {j q +i(0> Jq+iiv)} = J?([£j ??]), 77 G T. Introducing the Spencer operator 



D : R q+1 ->T* ® R q : (, q+1 ji(£ g ) - £ q+1 with local components {d^ k - £f , d l £^ - ...), we 



obtain the following general formula at order q: 

[€q,Vq] = {Zq+liVq+l} + i(0 D Vq+l ~ i{v)D£ q +l , V£ 9 , 77, £ R q 

where i() is the interior multiplication of a 1-form by a vector, and we let the reader check that 
such a definition no longer depends on the "lifts" t] q+ i over £, q7 r] q . 

EXAMPLE 3.3: (Affine transformations) n = 1, q = 2, X = R 3 

With evident notations, the system of finite Lie equations is defined by the single second order 
linear OD equation y xx = 0. Similarly, the sections of i?2 are defined by £ X x{%) = 0. Accordingly, 
the components of [£2,772] at order zero, one and two are defined by the totally unusual successive 
formulas: 



EXAMPLE 3.4: (Projective transformations) n = 1, q = 3, X = R 3 

The system of finite Lie equations is defined by the single third order nonlinear OD equation 
(yxxx/Ux) — ^(yxx/yx) 2 = and the sections of R3 are defined by £ xxx (x) = 0. The formulas for 
the bracket of Lie algebroid [R 3 , R 3 ] C R3 can be derived similarly but involve many more terms. 

EXAMPLE 3.5: (Volume preserving transformations) n arbitrary, q = 1,X = R n 

The only non-linear finite Lie equation is d(y 1 , ...,y n ) /d(x 1 , ...,x n ) — det(y^) — 1 and the sections 

of i?i are defined by the single relation Q — 0. Accordingly, we obtain: 



When summing on k and i, the first two terms disappear (as in Example 3.3 !) and we get therefore 
[Ri, Ri] C R\. We invite the reader to compare this result with the usual way on solutions where 
one defines 6 as the kernel of the Lie derivative of the (volume) n-form co — dx 1 A ... A dx n 

with respect to £ and then uses the well known formula [£(£), j0(rj)]uj = £([£, rj\)u), V£, 77 e T in 
order to obtain [O, 0] C O. 

Introducing a basis of AT* made by the dx 1 — dx 11 A...Adx lr with I = {i\ < ... < i r ), we may 
define the exterior derivative d : AT* — * A r+1 T* by setting u> — ujjdx 1 — > du; = diUJjdx 1 A dx 1 and 
one easily checks d 2 = d o d = 0. This is sufficient in order to define the (canonical linear) gauge 
sequence: 



It follows that £. 




([6,%]). 



{[£um])i - - VIZ! + CdrV* - ^drg 



AT* ® g A AT* ® Q A AT* ® 
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which is just the tensor product by Q of a part of the Poincare sequence for the exterior derivative. 

However, there are other differential sequences to be found in the literature and that we did not 
speak about, namely the Janet sequence, which is for sure the best known differential sequence, 
and the Spencer sequence. For short, if E,F, ... denote vector bundles over X, we use the same 
letters for the corresponding sets (sheaves to be exact) of sections and such an interpretation must 
be used whenever operators are involved. Starting from a vector bundle E (for example T) and a 
linear differential operator V : E — > F : £ — > 77, if we want to solve the linear system with second 
member X>£ = r\ even locally, one needs "compatibility conditions" (CC) in the form T>irj = 0. 
Denoting now F by F , we may therefore look for an operator T>i : F — ► F\ : 77 — > £ and so on. 
Under the assumption that V is involutivc (!), the french mathematician M. Janet has proved in 
1920 that such a chain of operators ends after n steps and we obtain the (canonical linear) Janet 
sequence, namely [8,22]: 

- 6 -> E 3. F S Fi 3 ... ^? F n - 

where, apart from the first operator which is of order q, the n operators following it are first order 
and involutive. The (canonical linear) Spencer sequence is the Janet sequence for the first order 
Spencer form, namely: 

where Co = R q and the first order involutive operators Di are induced by the Spencer operator. 
It follows that we only have at our disposal for any application where group theory seems to be 
involved, three linear differential sequences, namely the Janet sequence, the Spencer sequence and 
the gauge sequence. As these sequences are made by quite different operators, the use of one 
excludes the use of the others. 

In order to escape from this dilemna and for the sake of clarifying the key idea of the broth- 
ers Cosserat by using these new mathematical tools, we shall explain, in a way as elementary as 
possible while using only the linear framework, why the Janet sequence and the gauge sequence 
cannot be used in continuum mechanics. By this way we hope to convince the reader about the 
need to use another differential sequence, namely the SPENCER SEQUENCE, though striking it 
could be. Also we shall use very illuminating examples in order to illustrate our comments. 
First of all we exhibit the isomoprphism existing between the linear gauge sequence and the linear 
Spencer sequence. For this, if now G acts on X with a basis £ r = {^.dk} of infinitesimal generators, 
we may introduce the map: 

A°T* ® G -> J q (T) : X T (x) - X T (x)d^(x) 

It is known [21 ,p. 308] that this map becomes injective for q large enough and we may call R q its 
image for such a q. It follows from its definition that R q ~ R q +i is a system of infinitesimal Lie 
equations of finite type and we get for the Spencer operator (care to the notation) [22,26]: 

D : R q+1 T* <g> R q : £ q+1 - {8$ - - 8^ (x)8^(x)) 

Accordingly, the linear gauge sequence is isomorphic to the linear Spencer sequence: 

-> 9 -> A°T* ® R q 5. A 1 T* <g> R q Z A 2 T* ® R q 

the three isomorphisms being induced by the (local) isomorphism X x Q — > R q just described 
above. It is essential to notice that, though the linear Spencer sequence and the isomorphisms 
crucially depend on the action, by a kind of " miracle" the linear gauge sequence no longer depends 
on the action (look at the kernel of the first " d" ) . 

This result proves that the linear Spencer sequence generalizes the linear gauge sequence, with the 
major gain that it can be used even for Lie pseudogroups of transformations that are not coming 
from Lie groups of transformations. 

REMARK 3.6: When n = 3 and we deal with the Lie group of rigid motions, though surpris- 
ing it may look like, the formal adjoint of the (first) Spencer operator is exactly describing the 
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so-called stress and couple-stress equations found by the brothers Cosserat as we shall explain 
below. Though this result is in perfect agreement with the piczzoclectric or photoelastic coupling 
of elasticity and elcctromagnetism [23], it CONTRADICTS gauge theory where the lagrangians 
arc functions on A 2 T* <g> Q and NOT on T* ® Q as above [29]. 

Let us consider a (finite) volume J v dV in M 3 limited by a (closed) surface S — j s dS and let 
us introduce the outside unit normal (pseudo) vector n = (rij) on S. Let us now suppose that the 
surface element dS is acted on by the outside with a force dF = adS and a couple dC = jldS, 
where both a and /2 linearly depend on n through the stress tensor density a — (ct 4 - 7 ) and the 
couple-stress tensor density fi — (/i"- 7 = — ^ rjl ). It must be noticed that, using the standard 
Cauchy tetrahedral device, there is no reason " a-priori' to suppose that the stress tensor is sym- 
metric. We also suppose that the volume element dV is acted on by (see later on for the sign) a 
force —fdV and a momentum —rhdV with / = (f J ) and m — (m lJ — —m?' L ). 
Our purpose is now to study the equilibrium of the corresponding torsor fields with respect to an 
arbitrary cartesian frame Ox 1 x 2 x 3 . 

The equilibrium of forces is satisfied if we have the relation: 

f SdS- j fdV = =► / o^ndS - [ PdV = 
J s J V J s J V 

Using Stokes formula, this is equivalent to the well known stress equations: 

d l a^ = p 

This result shows that the surface density of forces a is equivalent, from the point of view of force 
equilibrium, to a volume density of forces / and this interpretation explains the sign adopted. 
Finally, the equilibrium of forces being satisfied, it is known that the equilibrium of momenta 
is also satisfied if it is satisfied with respect to an arbitrarily chosen cartesian frame. Hence, in- 
troducing the vector r— (x 1 , x 2 , x 3 ), the equilibrium of momenta is satisfied if we have the relation: 



/ (p+rAa)dS- [ (m + f A f)dV = 
J s J v 



s 

Projecting onto the axis Ox 3 , we obtain: 

/ (n r ' 12 + zV 2 - x 2 a rl )n r dS - [ (m 12 + x 1 f 2 - x 2 f x )dV = 
J s J v 

Using again Stokes formula and the previous stress equations, we obtain the couple-stress equations: 

d r n r ' ij + <r I] - a 31 = m 11 

This result shows that the surface density of forces a and couples ft is equivalent, from the point 
of view of torsor equilibrium, to a volume density of forces / and to a volume density of momenta 
to, provided the preceding stress and couple-stress equations are satisfied, and this interpretation 
explains the sign adopted. 

The combination of the stress and couple-stress equations have first been exhibited by E. and F. 
Cosserat in 1909 [4,5,pl37] without any static equilibrium experimental background and we now 
explain the key argument leading to the same equations just from group theoretical arguments. 
Of course, most of the engineering continua such as steel, concrete, glass, wate,... have the specific 
"constitutive laws" [i — 0, to = and we obtain therefore — a^ 1 , that is the stress tensor is 
symmetric, a situation not always encountered in liquid crystals. 

First of all, for the reader not familiar with the Spencer operator, we exhibit a similar result 
in a quite simpler 1-dimcnsional situation that will allow to recapitulate all the previous results.. 

EXAMPLE 3.7: Let us consider the Lie group of afunc transformations of the real line defined 
by the group action y = a x x + a 2 . The corresponding 2-dimensional Lie group G has coordinates 
a = (a 1 , a 2 ) and the group composition law is ab — (a 1 , a 2 )(6 1 , 6 2 ) = (a 1 ^ 1 ,^^ 2 + a 2 ) with inverse 
law a -1 = (1/a 1 , — a 2 /a 1 ). A basis of infinitesimal generators of the action may be obtained with 
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£ x = x-§^,£, 2 = ■§£ and we have [t- 1 , £ 2 ] = — £ 2 . As the exterior derivative is a linear involutive first 
order homogeneous operator, both with its formal adjoint, it could not have anything to do with 
the stress and couple-stress equations previously exhibited. 

Let us now deal with the Spencer sequence instead of the gauge sequence in this framework. 
First of all, we may consider the above Lie group of transformations as a Lie pseudogroup defined 
by the second order system of finite Lie equations y xx = 0. The corresponding system R 2 C J2{T) 
of infinitesimal Lie equations is £ xx — and the isomorphism between the gauge sequence and the 
Spencer sequence is induced by the map: 

(\ 1 (x),X 2 (x)) - (£(x) = xX\x) + X 2 (x),£ x (x) = X 1 (x),^ xx (x) = 0) 
The only two non-zero components of the Spencer operator become: 

d x £{x) -£x(x) = xd x X 1 (x) + d x X 2 (x),d x £ x (x) - = d x £_ x (x) = d x X 1 (x) 
Equating to zero these two components amounts to have: 

d x X 1 = , d x X 2 = 

Accordingly, gauging A just amounts to choose an arbitrary section of R 2 . 

The final touch, that could not be in the mind of any reader even on this very simple example, is to 
work out the formal adjoint of the Spencer operator. For this, multiplying the first component by 
a test function a(x), the second by a test function fi(x), then summing and integrating by parts, 
we get the following operator with second members (/, m) : 

d x a = f , d x [i + a = m 

The comparison with the previous mechanical results needs no comment. 

Taking into account this example, we now study the foundation of elasticity theory and we 
restrict the study to 2-dimcnsional (infinitesimal) elasticity for simplicity as the general situation 
has already been treated elsewhere and we just want to explain why the only founding problem of 
elasticity is the choice of an underlying Lie pseudogroup and an adapted differential sequence. 

1) The gauge sequence cannot be used: 

Looking at the book [5] written by E. and F. Cosserat, it seems at first sight that they just construct 
the first operator of the gauge sequence for n — 1 ([5], p 7), n — 2 ([5], p 66), n = 3 ([5], p 123) and 
finally n = 4 ([5], p 189) in the linearized framework. This is not true indeed because, according 
to the comment previously done, the adjoint operator is a divergence like operator, a situation not 
met in the couple-stress equations which is a linear operator with constant coefficients, but not 
of divergence type. In fact, a carefull study of the book proves that somewhere the action of the 
group on the space is used, but this is well hidden among many very technical formulas (Compare 
[5], p 136 with [21] p 295). 

2) The Janet sequence cannot be used: 

This result is even more striking because all texbooks of elasticity use it along the same scheme 
that we now describe. Indeed, after gauging the translation by defining the " displacement vector" 
£ = (^(x), £ 2 (x)) of the body, from the initial point x = (x 1 ^ 2 ) to the point y = x + £(x), 
one introduces the (small) "deformation tensor 3 ' e = l/2£(£)u as one half the Lie derivative with 
respect to £ of the cuclidcan metric oj, namely, in our case, the three components only (care): 

e = (en = e 12 = e 21 = 1/2(8^ + c^ 1 ), e 22 = d 2 £ 2 ) 

From the mathematical point of view, one uses to consider the Lie operator £>£ = : T —> S 2 T* 

(symmetric tensors), sometimes called Killing operator, through the formula: 

= LOrjdi£ r + LOi r dj^ r + £ r ' d r u)ij = Qjj = 2eij 
One may check at once the only generating "compatibility condition" Vie = 0, namely: 

<9ne 22 + <9 22 £n - 2<9i 2 ei 2 = 
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which is nothing else than the Riemann tensor of a metric, linearized at u>. 

However, the main experimental reason for introducing the first operator of this type of Janet 
sequence is the fact that the deformation is made from the displacement and first derivatives but 
must be invariant under any rigid motion. In the general case it must therefore have (n + n 2 ) — 
(n + n(n — l)/2) = n(n + l)/2 components, that is 3 when n — 2, and this is the reason why 
introducing the deformation tensor e. For most finite element computations, the action density 
(local free energy) w is a (in general quadratic) function of e and people use to define the stress 
by the formula cr y = dwjdtij which is not correct because w only depends on en, £12,622 when 
n = 2 as the deformation tensor is symmetric by construction. Finally, textbooks escape from this 
trouble by deciding that the stress should be symmetric and this is a vicious circle because we have 
proved it was not an assumption but an experimental result depending on specific constitutive 
laws. Accordingly, when n = 2, we should have tT y 'ey = cr 11 en + (2<r 12 )ei2 + cr 22 e22- Hence, even 
if we find the correct stress equations with this convenient duality keeping the factor " 2" , we have 
no way to get the stress and couple-stress equations together. 

3) Only the Spencer sequence can be used: 

Let us construct the formal adjoint of the Spencer operator by multiplying all the (2 x 2) + 2 = 6 
linearly independent nonzero components by corresponding test functions. For simplifying the 
summation, we shall raise and lower the indices by means of the (constant) euclidean metric, set- 
ting in particular ^ = u>i r £, r and £jj = uJi r £,j ■ The only nonzero first jets coming from the 2x2 
skewsymmetric infinitesimal rotation matrix of first jets are now £1,2 = —£2,1 while the second or- 
der jets are zero because isometries are linear transformations. We obtain in the present situation: 

<7 n di6 + a 12 (d^ 2 - £1,2) + o- 21 (d 2 ii - 6,1) + <r 22 &6 + M r,12 dr£i,2 

Integrating by parts and changing the sign, we just need to look at the coefficients of £1,^2 and 
$1,2, namely: 

6 — > fta 11 + d 2 a 21 = f 1 
6 — cV 12 +cV 22 = / 2 
$1,2 — > d r n r ' 12 + a 12 - a 21 = m 12 

in order to get the adjoint operator ad(D) : a" _1 T* ® R\ — » A n T* ® R\ : (a, fi) — > (/, to) relating 
for the first time the torsor framework to the Lie coalgebroid (see the beginning of section 3). 
These equations are exactly the three stress and couple-stress equations of 2-dimcnsional Cosserat 
elasticity. In the n-dimensional case, a similar calculation, left to the reader as an exercise of 
indices, should produce exactly the n(n + l)/2 stress and couple-stress equations in general. It is 
now possible to enlarge the group in order to get more equations, that is as many equations as the 
number of group parameters. Using the conformal group of space-time, the 4 elations give rise to 
4 nonzero second order jets only which allow to exhibit the 4 Maxwell equations for the induction 
(H, D) along lines only sketched by H. Weyl in [20] because the needed mathematics were not 
available before 1970. But, as we already said, this is another story !. 

REMARK 3.8: It becomes now clear that the dim(A 2 T* ® Q) = {n{n - l)/2)(n(n + l)/2) = 
n 2 (n 2 — l)/4 first order compatibility conditions for the Cosserat fields [3,4,5,11,27] (the so-called 
torsion and curvature of E. Cartan [2]) are described by the second (first order) Spencer operator 
in the Spencer sequence while the n 2 (n 2 — 1)/12 second order compatibility conditions for the de- 
formation tensor (the so-called Riemann curvature) are described by the second operator T>\ in the 
Janet sequence (see [22], Example 10, p 249-258). Accordingly, the torsion+ curvature of Cartan 
is not at all the generalization of the curvature of Riemann, contrary to what is still claimed in 
mathematical physics today. 

4) PARAMETRIZATION: 

The main tool in this last section will be duality theory, namely the systematic use of the formal 
adjoint of an operator. For this, if E is a vector bundle, we introduce its dual E* to be the vector 
bundle with inverse transition matrix (for example T* is the dual of T). The formal adjoint of an 
operator V : E -> F is the operator ad(V) : A n T* <x> F* -> A n T* ® E* defined by the relation: 
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< A, D£ >=< ad(V)\, £ > +duj 

where A is a test vector density and u> 6 /\ n ^ 1 T* comes from Stokes formula of integration by 
part. Any operator can be considered as the formal adjoint of another operator because we have 
the identity ad{ad{V)) = V. Also, if P, Q e D and P = a^d^, then ad(P) = X^-l)!^^ and 
ad(PQ) = ad(Q)ad{P). 

DEFINITION 4.1: We now define N from ad(D) exactly as we already defined M from V, that 
is N = DX/ Dad{T>)\. Of course N highly depends on the presentation of M. 

The following nontrivial theorem, first obtained in [18], provides a purely formal test for decid- 
ing about the existence of a parametrization and exhibiting one. At the same time it establishes a 
link with the two previous sections and is already implemented on the computer algebra package 
http://wwwb.math.rwth-aachen.de/OreModules 



THEOREM 4.2: A test for checking that a given differential module M is torsion- free proceeds 
in 5 steps: 

1) Write the defining system as the kernel of an operator V. 

2) Construct its formal adjoint ad(T>). 

3) Work out generating CC for ad(T>) as an operator ad(T)-i). 

4) Construct ad{ad{V_ 1 )) = D_ x . 

5) Work out generating CC for 2?_i as an operator V' . 

Then M is torsion- free if and only if V and T>' have the same solutions (both provide M). 

REMARK 4.3: We have ad(V oT>_ 1 ) = ad(X>_i) o ad(D) = and thus VoV-x = 0. Accordingly 
V ranges among the CC of 2?_i, a result symbolically written as V < V , and M is torsion- free if 
and only if V = V . In that case, with a slight abuse of language, the kernel of T> is parametrized by 
the image of X>_i. Otherwise, any CC in V but not in V provides a torsion element of M and £>_i 
provides a parametrization of the system determined by M/t(M) or, equivalently, M' — M/t(M) 
is the torsion- free module determined by T>' . 

For the reader aware of homological algebra, the following rather " magic!'' corollary depends on 
difficult technical though classical results of homological algebra ([1] but the best "basic" references 
are [7], Lemma 3.8, p 147; [16], Theorem 9, p 133; [25], Corollary 6.18, p 186). 

COROLLARY 4.4: The differential module extp(N, D) = t(M) does not depend on the presen- 
tation of N. 

EXAMPLE 4.5: As a first striking result, that does not seem to have been noticed by mechani- 
cians up till now, let us consider the situation of classical elasticity theory where V is the Killing 
operator for the euclidean metric, namely 2?£ = C(£)lo and T>\ the corresponding CC, namely the 
linearized Riemann curvature with n 2 (n 2 — 1)/12 components. In that case, as it is well known 
that the Poincare sequence for the exterior derivative is self-adjoint up to sign (for n = 3 the 
adjoints of grad, curl, div are respectively div, curl, grad) then the first extension module does not 
depend on the differential sequence used and therefore vanishes. Accordingly, ad{T>) generates the 
CC of ad{T>i). Hence, in order to parametrize the stress equations, that is ad(T>), one just needs 
to compute ad(T>i). For n — 2, we get: 

<K9ii£22 + ^22611 - 2<9i 2 ei 2 ) = <9 2 2</>eii ~ 2<9i 2 0ei2 + <9ii</>e 2 2 + <9i(—) + <9 2 (...) 

and recover ... the parametrization by means of the Airy function in a rather unexpected way. 
Exhibiting a parametrization for n > 3 thus becomes a straightforward exercise of computer alge- 
bra, the number of (pseudo)-potentials being n 2 (n 2 — 1)/12 along the last remark of the preceding- 
section. 

EXAMPLE 4.6: Contrary to the preceding situation, the second order Einstein operator V 
in General Relativity has n/(n + l)/2 components and is self-adjoint as it comes from Hubert 
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variational calculus of the scalar Ricmann density. The corresponding CC are the well known 
n divergence conditions and the adjoint 2?_i is therefore the Killing operator for the Minkowski 
metric. According to the previous example, its CC T>' is thus the linearized Riemann curvature 
with n 2 (n 2 — 1) /12 components. It follows that V < V strictly when n > 4 and no parametrization 
can be found in that case. 

EXAMPLE 4.7: We finally treat the case of Cosserat equations. For this, instead of using 
the Janet sequence as before, we now use the Spencer sequence which is isomorphic to the gauge 
sequence though with quite different operators. However, according to the general theorems of 
homological algebra, the existence of a parametrization does not depend on the differential se- 
quence used and therefore follows again, like in the first example of this section, from the fact 
that the Poincare sequence is self-adjoint up to the sign. In the present situation, we have 
C r = A r T* ® Ri ~ AT* ® Q with dim(Q) = n(n + l)/2. As we have shown in the last section 
that the Cosserat equations were just ad{D\), their first order parametrization is thus described 
by ad(D 2 ) and needs dim(C 2 ) = n 2 (n 2 — l)/4 (pseudo)-potentials according to the last remark of 
the preceding section. We provide the details when n — 2 but we know at once that we must use 
3 (pseudo)-potentials only. The case n = 3 could be treated similarly and is left to the reader as 
an exercise. 

The Spencer operator D\ is described by the equations: 

0i£i = An, c?i £ 2 -S,i,2 = A 12 ,d 2 £i -£ 2 ,i = A 2 i,d 2 & = A 22 ,0i£i,2 = B 1 ,d 2 £,i, 2 = B 2 

because R\ is defined by the equations £1^ = 0,£i i2 + £2,1 = 0,^2,2 = 0. 
Accordingly the 3 CC describing the Spencer operator D 2 are: 

02^11 - 01^21 + B 1 = 0, 2 Ai 2 - 0iA 22 + B 2 = 0, d 2 B x - d x B 2 = 

Multiplying these equations respectively by 1 , 2 , <f> 3 , then summing and integrating by part, we 
get ad(D 2 ) and the desired first order parametrization in the form: 

a 1 ' 1 = -020V 1 - 2 - -020V 2 ' 1 = 0i0V 2 > 2 - 0i0V M2 = -020 3 + 0V 2 ' 12 - 9i0 3 + 2 

as announced previously. It is important to notice that such a parametrization, which could also 
be obtined by localization, is coherent with the classical one already obtained by localization in 
Example 2.6 and which can be recovered if we cancel the couple-stress. 

5) CONCLUSION: 

We have already proved in our book "Lie Pseudogroups and Mechanics" (1988) that the group 
foundation of elasticity pioneered by E. and F. Cosserat (1909) was just described by adding to 
the non-linear Spencer sequence (1972) a convenient variational calculus. As a crucial conclu- 
sion and similarly to classical elasticity, even if the initial framework is non-linear, the resulting 
stress/couple-stress equations are linear. Accordingly, in order to parametrize these equations, 
one just needs to refer to infinitesimal Lie equations and the three corresponding canonical linear 
differential sequences that can be constructed, namely the gauge sequence, the Janet sequence and 
the Spencer sequence, though only the last one is useful. 

The main result of this paper, coming from unavoidable arguments of homological algebra, is that 
the Cosserat equations are just described by the formal adjoint of the first Spencer operator while 
their parametrization is provided by the formal adjoint of the second Spencer operator. This result 
finally explains why the parametrization of the stress/couple-strcss equations is first order while 
the parametrization of the classical stress is second order as it comes from dualizing the Janet 
sequence. 

We are happy to to pay such a tribute to E. and F. Cosserat on the occasion of this anniversary 
and hope it will give a new impulse for generalizing their work. 
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